clear sample;
for i=1:2
figure;
file = ['Data', num2str(i), '.csv']
file = readtable(file);
structarray = table2struct(file);
sample = [(0:1:255) ; zeros(1,256)]';
diff_var = [(0:1:255) ; zeros(1,256)]';
mobility = [(0:1:255) ; zeros(1,256)]';
diff_mobility = [(0:1:255) ; zeros(1,256)]';
formfactor = [(0:1:255) ; zeros(1,256)]';
freq = zeros(129,1);
for j = 0:3
j
[sample, freq(:,j+1)] = electrodeplot(j, structarray, sample);
electrode_mean(i,j+1) = mean(sample(:,j+1))
electrode_var(i,j+1) = std(sample(:,j+1))
electrode_diff = diff(sample(:,j+1))
diff_var(:,j+2) = std(electrode_diff)
mobility(:,j+2) = diff_var(:,j+2)/electrode_var(i,j+1)
diff_mobility(:,j+2) = [diff(mobility(:, j+2)) ; 0]
for f=1:256
formfactor (f,j+2) = (diff_mobility(f,j+2))/(mobility(f,j+2))
end
approxEnt = approximateEntropy(sample(:, j+1))
end
extractwave(1, 4, sample, "Delta");
extractwave(4, 7, sample, "Theta");
extractwave(8, 12, sample, "Alpha");
extractwave(12, 30, sample, "Beta");
extractwave(30, 100, sample, "Gamma");
end